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ABSTRACT 
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Interstellar polarization from aligned dust grains can be measured both in trans- 
P\| ■ mission at visible and near-infrared wavelengths and in emission at far-infrared and 

sub-mm wavelengths. These observations can help predict the behavior of foreground 
contamination of CMB polarimetry by dust in the Milky Way. Fractional polarization in 
emission from aligned dust grains will be at the higher range of currently observed val- 
ues of 4-10%. Away from the galactic plane, fluctuations in Q and U will be dominated 
by fluctuations in intensity, and less influenced by fluctuations in fractional polarization 
• and position angle. 
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1. Some Basics 
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Polarization due to magnetically aligned dust grains in the Milky Way galaxy is subject to both 
Q H i coherent and random processes (see Jones, 1996 for a review of grain alignment). Most observations 

of interstellar polarization have been at visual wavelengths, with a modest contribution at near- 
infrared (NIR) wavelengths. At these wavelengths the polarization is due to selective extinction by 
aspherical dust grains. In the far-infrared (FIR), these same grains are in emission. This emission 
is polarized and consequently a foreground contaminant for CMB polarimetry. A few simple facts 
about interstellar polarization should be noted. 

• Fractional polarization P is very probably NOT a function of magnetic field strength (Jones, 
Klebe & Dickey 1991). 

• Polarization in emission (P e ) behaves differently than polarization in transmission (P a )- 

• P a starts low, then increases with increasing optical depth and contains information on the 
random component of the aligning magnetic field. 

• P e starts high, then decreases with increasing optical depth. CMB polarimetry will sample 
the low optical depth regime. 



CMB polarization measurements will be concerned with contamination from fluctuations in 
Stokes vectors Q and U caused by variations in interstellar polarization in emission. Recall that Q 
and U are intensities, not fractions. Specifically: 
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Thus, Variations in Q and U can be due to fluctuations in fractional polarization, fluctuations in 
position angle, and fluctuations in total intensity. For example, a region of sky with the same frac- 
tional polarization in emission (uniform grain alignment across the sky) will still show fluctuations 
in Q and U if the intensity is varying across the region. At high latitudes, this effect will turn out 
to be the major source of contamination. 

Polarimetry at FIR and sub-mm wavelengths would seem the best source of observations 
for characterizing the CMB polarized foreground. These observations measure polarized emission 
from grains, and in the sub-mm extend to wavelengths close to those to be observed in upcoming 
and proposed CMB polarimetry experiments. However, FIR and sub-mm observations require 
significant optical depth to achieve detectable emission, and consequently have largely been confined 
to dense regions in the galactic plane. At high latitudes, where CMB polarimetry will be most 
effective, the emission is very optically thin and nearly impossible to observe in polarization. The 
reader should see Hildebrand et al. (2000) for a thorough discussion of FIR polarimetry. 



Optical polarimetry can, however, measure polarization in transmission (P a ) at high latitudes 
because dust extinction is so strongly wavelength dependent. We need to relate optical, NIR 
and FIR polarimetric observations in order to predict the characteristics of foreground polarized 
emission at high latitudes. This work is in its infancy, but some useful conclusions can be made. 

Jones, Klebe &; Dickey (1991) present a simple model of grain alignment by the interstellar 
magnetic field that they use to explain the general trend of polarization (P a ) with optical depth 
(t) observed at optical and NIR wavelengths. They assume the magnetic field is strong enough 
everywhere to fully align the grains and that any departure from a simple linear relation between 
P a and r is due to fluctuations in the magnetic field geometry. Their results are shown in Figure 
1. The upper model line is the case for a perfectly uniform magnetic field geometry, the lower 
model line is for a totally random magnetic field geometry. The data clearly lie between these 
two extremes. The solid middle line is the model result with a 50/50 mix of random and uniform 
components and a decorrelation length for the random component of tk = 0.1. 



2. 



Observations of Interstellar Polarization 
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Fig. 1. — Fractional polarization vs. optical depth at K (2.2 fim) and the JKD model. The solid model line 
corresponds to a 50/50 mix of random and constant components to the interstellar magnetic field with a decorrelation 
length of tk = 0.1 

Similar results have been obtained by Fosalba et al. (2002) using a much larger optical data 
set. These model results show that fluctuations in P a will be strongly influenced by fluctuations in 
the magnetic field geometry only for optical depths greater than r ~ 0.03 — 0.1 at 2.2 /im (about 
Ay = 0.3 — 1.0). At high galactic latitudes, where the optical depth is much less than this, we 
would expect the magnetic geometry to be relatively smooth, and therefore not a significant factor 
influencing fluctuations in Q and U in emission at FIR wavelengths. 

This conclusion can be seen more clearly by using the parameters from the JDK model to 
predict the polarization in emission. The results of this exercise are shown in Figure 2, where the 
predicted fractional polarization at FIR wavelengths is plotted vs. the number of decorrelation cells 
along the line of sight. The model computes the predicted average polarization for an ensemble of 
lines of sight. Note that P e is highest at the low optical depths, since all of the grains are emitting 
with the same alignment geometry. This is the regime in which CMB observations will be most 
useful. 

Based on the results from NIR polarimetry (Figure 1), one decorrelation cell corresponds 
roughly to Ay = 0.3 — 1, or in physical distance units about 200-500 pc in the diffuse ISM. 
This corresponds to a greater optical depth than seen at high latitudes and one such cell would 
cover several degrees across the sky. Thus, the JKD model predicts that the CMB foreground 
contamination will most likely have a high intrinsic fractional polarization but a relatively uniform 
projected geometry over many degrees at high latitudes. 

In Figure 2 we have assumed the maximum polarization of completely aligned grains is 7%, a 
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Far Infrared Polarized Emission, Optically Thin Case 
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Fig. 2. — Predicted average fractional polarization in the FIR vs. number of decorrelation cells using the JKD 
model. The middle line corresponds to a 50/50 mix of random and constant components. The polarization for this 
model is less than the maximum (7%) due to the fact that the random component of the magnetic field can point in 
any direction, including along the line of sight. 

tyical upper bound for FIR measurements in dense star forming regions. More recent observations 
in regions with lower optical depth and observations at sub-mm wavelengths strongly suggest that 
10% would be a more appropriate choice for the maximum fractional polarization warm dust would 
have at the wavelengths of CMB experiments (Hildebrand et al. 2000). 

Direct optical polarimetry in regions of the sky a few degrees across at high latitudes confirms 
this prediction. For example, the region from I = 275° — 295°, b = 60° — 70° shown in Figure 3 has 
optical polarimetry for over 25 stars in the compilation by Heiles (1996). In this region we find a 
very uniform position angle (0 rms = 0.25rad). The observed fluctuations in fractional polarization 
(Prms/ < P >= 0.4) are about a factor of two higher than expected given the observed d rms . This 
is presumably due to fluctuations in optical depth across the region, not fluctuations in magnetic 
field geometry. 



This leaves fluctuations in total intensity as the major source of fluctuations in Q and U from 
emission of aligned dust at high latitudes. Using the 100 /xm map from Schlegel, Finkbeiner & Davis 
(1998) and the predicted average fractional polarization from Figure 2 for one decorrelation cell, we 
can predict the fluctuations in Q and U due to fluctuations in total intensity at FIR wavelengths. 



3. FIR Surface Brightness 
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Optical Polarimetry 
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Fig. 3. — Optical polarimetry of stars at high galactic latitude in a region taken from the compilation by Heiles 
(1996). The measured dispersions in position angle and fractional polarization are 9 rma — 0.25rad and P rmB j < 
P >= 0.4. Note how uniform the magnetic field geometry is over a span of several degrees. 

To make a prediction at a frequency of 150 GHz, we use the following input parameters: 

1. Assume < P e > and 9 are both constant with < P e >= 5%. 

2. The measured fluctuations of the intensity in the central 40° x 40° region of the 100 fim map 
is I rm s = 0.4MJy/sr using 10' bins. 

3. Extrapolate the intensity from 100 /im to 150 GHz using a power law index of 1.75 for the 
emissivity. 

4. Arbitrarily assign Q to lie along the magnetic field direction. In this configuration the average 
value of U will be 0. 

With these parameters, we find that Q rrns ~ 10 _ Icmb due to fluctuations in total intensity alone. 

We can improve this estimate somewhat if we include the lesser effects of fluctuations in the 
magnetic field geometry. Taking the observed dispersion in position angle for the region of sky 
shown in Figure 3 {9 rms = 0.25rad) and using the JKD model to predict the fluctuations in P a , 
we find P rms / < P >= 0.2, about half the observed value. Combining the effects of fluctuations in 
intensity at 150 GHz with fluctuations in magnetic field geometry for this one piece of the sky, we 
obtain the predicted values for Q rms and U rms given in Table 1. 
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Table 1. Predicted Contamination 



GHz 


Qrms/ ICMB 


Urms/ICMB 
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100 


0.25 


0.12 


150 


0.91 


0.44 


200 


2.6 


1.3 
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